Local structure study of In x Gai x As semiconductor alloys using High Energy 

Synchrotron X-ray Diffraction 
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Nearest and higher neighbor distances as well as bond length distributions (static and thermal) 
of the In x Gai_ x As (0 < x < 1) semiconductor alloys have been obtained from high real-space 
resolution atomic pair distribution functions (PDFs) . Using this structural information, we modeled 
the local atomic displacements in In x Gai_ x As alloys. From a supercell model based on the Kirkwood 
potential, we obtained 3-D As and (In,Ga) ensemble averaged probability distributions. This clearly 
shows that As atom displacements are highly directional and can be represented as a combination 
of (f00) and (111) displacements. Examination of the Kirkwood model indicates that the standard 
deviation (a) of the static disorder on the (In,Ga) sublattice is around 60% of the value on the 
As sublattice and the (In,Ga) atomic displacements are much more isotropic than those on the As 
sublattice. The single crystal diffuse scattering calculated from the Kirkwood model shows that 
atomic displacements are most strongly correlated along (110) directions. 
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I. INTRODUCTION 

Semiconductor alloys are known as technologically im- 
portant materials for their wide applications ia optoelec- 
tronic devices such as lasers and detectors.^ The local 
structure information is of fundamental importance in 
understanding the alloy systems because their physical 
properties are strongly influenced by the local atomic 
displacements present in the alloys. For example, it is 
known that chemical and compositional disorder strongly 
affect ,-the electronic structure of zinc-blende type al- 
loysB lI and their enthalpies of formation.ffl 

In this paper we present a detailed study of the lo- 
cal and average structure of the In x Gai_ x As alloy se- 
ries. The average structure of In x Gai_ x As was studied 
by Woolley.ty The structure is of the zinc blende type 
(F43m)EJ over the entire alloy range. The lattice param- 
eters, and therefore the average In-As and Ga-As bond 
lengths, interpolate linearly between the values of the 
end-members according to Vegard's law £3 However, con- 
sideration of the local structure reveals a very different 
situation. The local structure of In x Gai_ x As was first 
studied by Mikkelson and Boyce using extended x-ray 
absorption fine structure (XAFS).E3 According to this 
experiment, the individual nearest neighbor (NN) Ga- 
As and In-As distances in the alloys are rather closer to 
the pure Ga-As and In-As distances. Further XAFS ex- 
periments showed that this is quite geaepal behavior for 
many zinc-blende type alloy systems. Eilila Since then a 



number of theoretical and model studies have been car- 
ried out on the semiconductor alloys to understand—how 
the alloys accommodate the local displacementsJiZfO 

Until now these models and theoretical predictions are 
tested mainly by the comparison with XAFS data. The 
XAFS results give information about the nearest neigh- 
bor and next nearest neighbor distances in the alloys 
but imprecise information about bond-length distribu- 
tions and no information about higher-neighbor shells. 
This limited structural data makes it difficult to differ- 
entiate between competing models for the local struct 
ture. For example, even a simple radial force modeled 
rather accurately predicts the nearest neighbor distances 
of In x Gai_ x As alloys in the dilute limit. Therefore, 
one needs more complete structural information includ- 
ing nearest neighbor, far-neighbor distances, and bond 
length distributions to prove the adequacy of model struc- 
tures for these alloys. 

The atomic pair distribution function (PDF), G(r), 
measures the probability of-4inding an atom at a dis- 
tance r from another atom.t3 One of the advantages of 
the PDF method over other local probes such as XAFS 
is that it gives both local and intermediate range infor- 
mation because both Bragg peaks and diffuse scattering 
are used in the analysis. It is also possible to obtain in- 
formation about the static bond length distribution from 
the PDF peak width and about correlations of atom dis- 
placements.^ 

In this paper, we present a detailed X-ray diffraction 
study of In x Gai_ x As, (0 < x < 1). A preliminary analy- 
sis of the data has been published elsewhere.EJ Using high 
energy synchrotron x-rays, we measured the total scat- 
tering structure function, S(Q), of the In x Gai_ x As alloy 
system extended to high Q (Q m ax = 45 A -1 ) where Q is 
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the magnitude of the momentum transfer of the scattered 
x-rays (Q = 47rsin#/A for elastic scattering). From these 
structure functions we obtained the corresponding high 
real-space resolution PDFs through a Fourier transform 
according to 

2 I" 00 

G(r) = - / Q[S(Q) -1] sin QrdQ. (1) 
* Jo 

In these PDFs, the first peak is clearly resolved into two 
sub-peaks corresponding to the Ga-As and In-As bond 
lengths.LZl The evolution of the bond-length with doping 
gives good agreement with XAFS. For the far- neighbor 
peaks, the peak widths are much broader in the alloy 
samples than those of the pure end-members reflecting 
the increased disorder. We model the local structure of 
In x Gai_ x As alloys ifsjng a supercell modeled based on the 
Kirkwood potential^ which gives good agreement with 
the alloy data with no adjustable parameters. The results 
of the modeling have been analyzed to reveal the average 
atomic static distribution on the As and (In,Ga) sublat- 
tices. Finally, we have calculated the diffuse scattering 
that one would get from the Kirkwood model. This com- 
pares qualitatively well with-published diffuse scattering 
results from Ino.53Gao.47As.E3 

II. EXPERIMENTAL DETAILS 
A. Data collection 

The alloy samples, with compositions In x Gai_ x As, (x 
= 0, 0.17, 0.33, 0.5, 0.83, 1) were prepared by a melt and 
quench method. An appropriate fraction of InAs and 
GaAs crystals were powdered, mixed and sealed under 
vacuum in quartz ampoules. The samples were phjeated 
beyond the liquidus curve of the respective alloytlEl to 
melt them and held in the molten state for 3 hours before 
quenching them in cold water. The alloys were powdered, 
resealed in vacuum, and annealed just below the solidus 
temperature for 72-96 hours to increase the homogeneity 
of the samples. This was repeated until the homogeneity 
of the samples, as tested by x-ray diffraction, was satis- 
factory. X-ray diffraction patterns from all the samples 
showed single, sharp diffraction peaks at the positions 
expected for the nominal alloy similar to the results ob- 
tained by Mikkelson and Boyce.td 

High energy x-ray powder diffraction measurements 
were conducted at the A2 wiggler beamline at Cor- 
nell High Energy Synchrotron Source (CHESS) using in- 
tense x-rays of 60 KeV (A = 0.206 A). The incident x- 
ray energy was selected using a Si(lll) double-bounce 
monochromator. All measurements were carried out in 
flat plate symmetric transmission geometry. In order 
to minimize thermal atomic motion in the samples, and 
hence increase the sensitivity to static displacements of 
atoms, the samples were cooled down to 10 K using a 
closed cycle helium refrigerator mounted on the Hubcr 



6 circle diffractormeter. The samples were uniform flat 
plates of loosely packed fine powder suspended between 
thin foils of kapton tape. The sample thicknesses were 
adjusted to achieve sample absorption fj,t ~ 1 for the 
60 KeV x-rays, where /i is the linear absorption coeffi- 
cient of the sample and t is the sample thickness. 

The experimental data were collected up to Q m ax — 45 
A -1 with constant AQ steps of 0.02 A -1 . This is a very 
high momentum transfer for x-ray diffraction measure- 
ments. For comparison, Qmax from Cu K a x-ray tube is 
less than 8 A -1 . This high Q m ax is crucial to resolve the 
small difference (as 0.14 A) in the In-As and Ga-As bond 
lengths. 

To minimize the measuring time, the data were col- 
lected in two parts, one in the low Q region from 1 to 13 
A" 1 and the other in mid- high Q region from 12 to 50 
A -1 . Because of the intense scattering from the Bragg 
peaks, in the low Q region the incident beam had to be 
attenuated using lead tape to avoid detector saturation. 
The maximum intensity was scaled so that the count rate 
across the whole detector energy range in the Ge detec- 
tor did not exceed ~ 2 x 10 4 s" 1 . At these count-rates 
detector dead-time effects are significant but can be reli- 
ably corrected as we describe below. To reduce the ran- 
dom noise level below 1%, we repeated runs until the 
total elastic scattering counts become larger than 10,000 
counts at each value of Q. Also, to obtain a better pow- 
der average the sample was rocked with an amplitude of 
±0.5° at each Q-position. The scattered x-rays were de- 
tected using an intrinsic Ge solid state detector. The sig- 
nal from the Ge detector was processed in two ways. The 
signal was fed to a multi-channel analyzer (MCA) so that 
a complete energy spectrum was recorded at each data- 
point. The signals from the elastic and Compton scat- 
tered radiation could then be separated using software 
after the measurement. In parallel, the data were also 
fed through single-channel pulse-height analyzers (SCA) 
which were preset to collect the elastic scattering, Comp- 
ton scattering, and a wider energy window to collect 
both the elastic and Compton signals. For normaliza- 
tion, the incident x-ray intensity was monitored using an 
ion chamber detector containing flowing Ar gas. 

For the SCAs, the proper energy channel setting for 
the elastic scattering is crucial. Any error in the channel 
setting could cause an unknown contamination by Comp- 
ton scattering and make data corrections very difficult. 
There's no such problem in the MCA method since the 
entire energy spectrum of the scattered radiation is mea- 
sured at each value of Q. The main disadvantage of the 
MCA method is that it has a larger dead-time, although 
this can be reliably corrected as we show below. Fig. |l| 
shows a representative MCA spectrum taken from the 
InAs sample at Q — 45 A. It is clear that the Compton 
and elastic scattering are well resolved at this high mo- 
mentum transfer. The elastically scattered signal, which 
contains the structural information, is obtained by inte- 
grating the area under the elastic scattering peak. 
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FIG. 1. MCA spectrum of InAs at Q=45 A -1 . Peaks 
in the spectrum from the elastic and Compton scattering 
are labelled, as is a peak from an electronic pulser used for 
dead-time correction. The other peaks in the spectrum come 
from various fluorescence and escape peaks. 
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FIG. 2. Calculated double scattering ratio, I2/I1 where 
I2 is the intensity due to double scattering events and Ij the 
single scattering intensity, for the Ino.5Gao.5As alloy in trans- 
mission geometry with a fit = 1.32 appropriate for our sample. 



B. Data Analysis 

The pjeasured x-ray diffraction intensity may be ex- 
pressed^ by 

jmea [Q) = PA[N{I coh + jinc + jnrul }] ; (2) 

where P is the polarization factor, A the absorption fac- 
tor, N the normalization constant, and I^° h , P e " c , I™ ul 
are the coherent single scattering, incoherent (Comp- 
ton), and multiple scattering intensities, respectively, per 
atom, in electron units. The total scattering structure 
function, S(Q), is then defined as 

S(Q) = [I^ -((f) -iff)}/ iff, (3) 

where (/) = (f(Q)) is the sample average atomic form 
factor and (f 2 ) is the sample average of the square of 
the atomic form factor. Therefore, to obtain S(Q) from 
the measured diffraction data, we have to apply correc- 
tions such as multiple scattering, polarization, absorp- 
tion, Compton scattering and Laue diffuse corrections 
on the raw data-ErLJ 

The corrections were carriecLput using a home- written 
computer program, PDFgetXo that is able to utilize the 
MCA data. The results obtained using the MCA ap- 
proach are_very similar to those obtained using the SCA 
approach.EZl It appears that both approaches work well 
for quantitative high energy x-ray powder diffraction. 
One possible advantage of the MCA method is that en- 
ergy windows of interest can be set after the experiment 
is over which is precluded if data are only collected using 
SCA's. 

We briefly describe some of the features of the data 
correction using PDFgetX. Data are first corrected for 
detector dead-time. In this experiment, we used the 
pulser method.E3 A pulse-train from an electronic pulser 
of known frequency is fed into the detector preamp. The 
voltage of the pulser pulses are set so that the signal 
appears in a quiet region of the MCA spectrum. The 



measured counts in the pulser signal in the MCA (or, in- 
deed, in an SCA window set on the pulser signal) is then 
recorded for each data point. The data dead-time cor- 
rection is then obtained by scaling the raw data by the 
ratio of the known pulser frequency and the measured 
pulser counts. This method accounts for dead-time in the 
preamp, amplifier and MCA/SCA electronics but not in 
the detector itself. However, in general the dead-time is 
dominated by the pulse-shaping time in the amplifier or 
the analogue-digital conversion in the MCA or SCA and 
so this method gives rather accurate dynamic measure- 
ment of the detector dead-time. An alternative dead- 
time correction protocol for correcting MCA data is to 
use the MCA real-time/live-time ratio. This works rea- 
sonably well if the MCA conversion time is the dominant 
contribution to the detection dead-time. This approach 
gave similar results to the pulser correction in this case. 

Multiple (mainly double) scattering can be a problem 
if samples are relatively thick and the radiation is highly 
penetrating as in the present case. The multiple scat- 
tering contribution contains no usable structural infor- 
mation and must be removed from the measured inten- 
sity. It depends on sample thickness and many other 
sample dependent factors such as attenuation coefficient, 
atomic number and weight of sample constituent E3~t3 It 
increases as the sample becomes thicker in both trans- 
mission and reflection geometry. The multiple scatter- 
ing correction was-calculated using the approach sug- 
gested by WarrenESrEa in the isotropic approximation. 
Calculation of the multiple scattering intensity is consid- 
erably simplified when the elastic and Compton signals 
are separated as is done here since only completely elas- 
tic multiple scattering events need to be considered. In 
In x Gai_ x As samples, the multiple scattering ratio was 
around 10% maximum at high Q in transmission geome- 
try. This result suggests that the proper multiple scatter- 
ing correction becomes important in the high Q region. 
Fig. H shows the double scattering ratio calculated for 
the Ino.5Gao.5As sample. 

The x-ray polarization correction is almost negligi- 
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FIG. 3. Comparison between Compton and elastic scat- 
tering intensities measured in Ino.5Gao.5As. 



ble for synchrotron x-ray radiation because the incident 
beam is almost completely plane polarized perpendicu- 
lar to the scattering plane. As a result there is virtually 
no angle dependence, to the measured intensity due to 
polarization effects.EB 

The Compton scattering correction is very important 
in high-energy x-ray diffraction data analysis. It can be- 
come larger than the coherent scattering intensity at high 
Q, as is evident in Fig. |l|. Even a small error in deter- 
mining the Compton correction can lead to a significant 
error in the coherent scattering intensity in the high Q 
region. However, in this region of the diffraction pattern 
the elastic and Compton-shifted scattering are well sep- 
arated in energy and can be reliably separated using the 
energy resolved detection scheme we used here. At low-Q 
the Compton-shift is small and the Compton and elas- 
tic signals cannot be explicitly separated unless a higher 
energy resolution measurement is made, for example, us- 
ing an analyzer crystal. However, the Compton inten- 
sity is much lower and the coherent scattering intensity 
is much larger. In this region a theoretically calculated 
Compton signal can be subtracted from the data contain- 
ing both elastic and Compton scattering. Uncertainties 
in this process have a very small effect on the resulting 
S(Q). Fig. H shows the signals from the Compton and 
elastic scattering in the Ino.5Gao.5As sample. At low Q, 
some contamination from the elastic scattering is appar- 
ent in the Compton channel. For the Compton scattering 
correction, we followed two steps. In the high Q region 
the Compton scattered signal was directly removed by 
integrating a narrow region of interest in the MCA spec- 
trum which only contained the elastic peak. In the low 
Q region-.we calculated the theoretical Compton scat- 
tcringEjO and subtracted this from the combined (unre- 
solved) Compton plus elastic scattering signal. These two 
regions were smoothly interpolated using a window func- 
tion, following the method of Rulandcll in which the the- 
oretical Compton intensity is smoothly attenuated with 
increasing Q. 

At very high Q values, due to the Debye- Waller factor, 
the Bragg peaks in the elastic scattering signal disap- 
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FIG. 4. The reduced total scattering structure function 
[S(Q) — 1]Q for In^Gai-^As measured at 10K. The data-sets 
are offset for clarity. The high-Q region is shown on an ex- 
panded scale (x 3) to highlight the presence of diffuse scat- 
tering. 



pear and the normalized intensity asymptotes to (f(Q) 2 )- 
This fact allows us to obtain an absolute data normal- 
ization by scaling the data to line up with (f(Q) 2 ) in the 
high-Q region of the diffraction pattern. Finally, the to- 
tal scattering structure function, S(Q), is then obtained 
using Eq. || and the corresponding PDFs, G(r) are ob- 
tained according to Eq. |]. 



III. RESULTS 

Fig. |J shows the experimental reduced total scatter- 
ing structure functions, F(Q) — Q[S(Q) — I], for the 
In x Gai_ x As alloys measured at 10 K. It is clear that the 
Bragg peaks are persistent up to Q ~ 35 A -1 in the 
end-members, GaAs and InAs. This reflects both the 
long range order of the crystalline samples and the small 
amount of positional disorder (dynamic or static) on the 
atomic scale. In the alloy samples, however, the Bragg 
peaks disappear at much lower Q-values but still many 
sharp Bragg peaks are present in the mid-low Q region. 
Instead, oscillating diffuse scattering which contains local 
structural information is evident in high Q region. The 
observation of Bragg peaks reflects the presence of long- 
range crystalline order in these alloys. The fact that the 
Bragg peak intensity disappears at lower Q-values in the 
alloys than the end-members reflects that there is signifi- 
cant atomic scale disorder in the alloys as expected. The 
oscillating diffuse scattering in the high-Q region orig- 
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FIG. 5. The reduced PDF, G(r) for In^Gai-^As mea- 
sured at 10 K. The data-sets are offset for clarity. 

inates from the stiff nearest-neighbor In-As and Ga-As 
covalent bonds. 

Fig. H shows the corresponding reduced PDFs, G(r), 
obtained using Eq. |l|. In the alloys, it is clear that the 
first peak is split into a doublet corresponding to shorter 
Ga-As and longer In-As bonds.c3 The position in r of 
the left and right peaks does not disperse significantly 
on traversing the alloy series. This shows that the lo- 
cal bond lengths stay close to their end-member values 
and do not follow VegarrUs, law, in agreement with the 
earlier XAFSt2l and PDFc3 reports. However, already 
by 10 A the structure is behaving much more like the 
average structure. For example, the doublet of PDF 
peaks around 11 A in GaAs (Fig. |J) remains a doublet (it 
doesn't become a quadruplet in the alloys) and disperses 
smoothly across the alloy series to its position at around 
12 A in the pure InAs. This shows that already by 10 A 
the structure is exhibiting Vegard's law type behavior. 

It is also notable that for the nearest neighbor PDF 
peak, the peak widths are almost the same in both alloys 
and end-members but for the higher neighbors, the peaks 
are much broader in the alloys than in the end-members. 

IV. ATOMIC DISPLACEMENTS IN THE ALLOYS 
A. Modeling 

A simplified view of the structural disorder in the 
A x Bi_ x C type tetrahedral alloys can be intuitively vi- 
sualized by considering simple tetrahedral clusters cen- 
tered about C sites (the unalloyed site). In the random 




FIG. 6. Schematic diagram of As displacements in cluster 
(a) type II, (b) type III, and (c) type IV. Cluster types are 
discussed in the text. At the corner, large dark circle and 
small grey circle show In and Ga atom position respectively. 
At the center, the grey and dark circles correspond to the As 
atom position before and after displacement, respectively. 

alloy this site can have 4 A-neighbors (type-I), 3 A- and 
1 B-neighbors (II), 2 A- and 2 B- neighbors (III), 1 A 
and 3 B-neighbors (IV) or 4 B neighbors (V). We as- 
sume that the mixed site (A,B) atoms stay on their ideal 
crystallographic positions. By considering each cluster 
type in turn we can predict the qualitative nature of the 
atomic displacements present in the alloy. Let the A 
atoms be larger than the B atoms. In clusters of type 
I and V the C atom will not be displaced away from 
the center of the tetrahedron. As is shown in Fig. |^, in 
type II clusters the C atom will displace away from the 
center directly towards the B atom. This is a displace- 
ment in a (111) crystallographic direction. In type III 
clusters it will displace in a direction between the two B 
atoms along a (100) crystallographic direction. Finally, 
in type IV clusters it will again be a (111) type displace- 
ment but this time in a direction directly away from the 
neighboring A atom. Such a cluster model was used to 
make quantitative comparisons with the nearest neighbor 
bond distances observed in XAFS measurementslij over 
the whole alloy series. t3 However, it was later shown that 
the agreement was largely accidental due to the boundary 
conditions chosen that all the clusters should have the av- 
erage size determined from Vegard's law .113 Nonetheless, 
it is interesting to compare the prediction of this simple 
cluster model with the nearest-neighbor PDF peaks mea- 
sured here since, for the first time, we have an accurate 
measurement of the bond length distributions as well as 
the bond lengths themselves. 

Each cluster is independenthz_relaxed according to the 
prescription of Balzarotti et alt3 to get the bond- lengths 
within each cluster type. Assuming a random alloy the 
number of each type of cluster that is present can be esti- 
mated using a binomial distribution. This gives the static 
distribution of bond lengths predicted by the model. 
These are then convoluted with the broadening expected 
due to thermal motion. This was determined by mea- 
suring the width of the nearest neighbor peaks in the 
end-member compounds, InAs and GaAs. The result is 
shown in Fig. 0(a). It is clear that, although the cluster 
model gets the peak positions reasonably correct as epu 
cmplificd by the agreement it gets with XAFS data,Efl 
it does rather a poor job of explaining the shape of the 
measured pair distribution. The major discrepancy is 
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FIG. 7. Comparison between Experimental PDF (open 
circles) and model PDF (solid line) for Ino.5Gao.5As. (a) 
Tetrahedral cluster model with no disorder present on (In,Ga) 
sublattice. The sub-peaks represent the contributions from 
each type of cluster. Type I (x), type II (□), type III (O), 
type IV (A), and type V(*). (b) The model PDF is calculated 
in the Pauling limit. The peak positions were obtained from 
the InAs and GaAs bond lengths in the end-members (solid 
line) and the InAs and GaAs bond lengths in the Ino.5Gao.5As 
PDF (dashed line). See the text for details, (c) Kirkwood su- 
percell model. 



that too much intensity resides at, or close to, the undis- 
placed position leading to an unresolved broad first PDF 
peak in sharp contrast to the measurement. In contrast, 
we show in Fig. 0(b) the nearest-neighbor atomic pair 
distribution in the Pauling limit ,0 again broadened by 
thermal motion. The peak positions were obtained by 
using the bond- lengths of the end- member compounds. 
It is clear that this actually does a better job than the 
cluster model, though it slightly, and not surprisingly, 
overemphasizes the splitting. The dashed line in this fig- 
ure shows the peak profile that we obtain if we make 
the assumption that the nearest neightu 



jond length 

changes in the alloy as seen in the Z-plot,t3E3 but there is 
no increase in the bond length distribution. Again, this 
gives rather good agreement emphasizing the fact that 
there is very little inhomogeneous. strain to the covalent 
bond length due to the alloying. E3 nn 
A better model for the structure of these alloysO'O 
is obtained from a relaxed supejxell of the alloy sys- 
tem using a Kirkwood potential.^ The potential con- 
tains nearest neighbor bond stretching force constants 




FIG. 8. Comparison between Experimental PDF (open 
circles) and model PDF (solid line) for In x Gai_ x As. The 
model was the Kirkwood supercell model. The parameters a 
and j3 are refined from the end-members and the PDFs for 
the alloys shown here are then calculated with no adjustable 
parameters. 



a and force constants (3 that couple to the change 
in the angle between adjacent nearest neighbor bonds. 
In this relaxed supercell model, th£-. force constants 
were adjusted to fit the end-membersEll with aca-As = 

96N/m, ai n _A S = 97N/m, /3 G a-As-Ga = /^As-Ga-As = 

lON/m and An-A s -in = /3As-in-As = 6N/m. The ad- 
ditional angular force constants required in the alloy arc 
taken to be the geometrical mean, so that /3ca-As-in = 
v/^Ga-As-Ga-An-As-in)- The PDFs for the alloys could 
then be calculated in a self-consistent way for all the al- 
loys with no adjustable parameters. a In this model, the 
lattice dynamics are also included in a completely self- 
consistent way. Starting with the force constants and 
the Kirkwood potential, the thermal broadening of the 
PDF peaks at any temperature can be determined di- 
rectly from the dynamical matrix and this is how the 
PDFs were calculated in the present case£3 The model- 
PDF is plotted with the data in Fig. || with the nearest- 
neighbor peak shown on an expanded scale in Fig. 0(c). 
The excellent agreement with the data over the entire 
alloy range suggests that the simple Kirkwood potential 
provides an adequate starting point for calculating dis- 
torted alloy structures in these III-V alloys. Note that 
in comparing with experiment, the theoretical PDF has 
been convoluted with a Sine function to incorporate the 
truncation of the experimental data at Q ma x — 45 A. 
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FIG. 9. Iso-probability surface for the ensemble averaged 
As atom distribution. The surfaces plotted all enclose the 
volume where As atoms will be found with 68 % probability, 
(a) Ino.17Gao.83As (b) Ino.33Gao.67As (c) Ino.5oGao.5oAs (d) 
Ino.83Gao.17As. In each case, the probability distribution is 
viewed down the [001] axis. 



B. 3-D atomic probability distribution 

Now, we analyze the relaxed supercell of alloy system 
obtained using a Kirkwood potential to get the average 
three dimensional atomic probability distribution of As 
and (In,Ga) atoms. Fig. || shows iso-probability surfaces 
for the As site in the In x Gai_ x As alloy. The probability 
distributions were created by translating atomic positions 
of the displaced arsenic atoms in the supercell (20x20x20 
cubic cell) into a single unit cell. To improve statistics, 
this was done 70 times. The surfaces shown enclose a 
volume where the As atom will be found with 68 % prob- 
ability. The probability distribution is viewed down the 
[001] axis. It is clear that the As atom displacements, 
though highly symmetric, are far from being isotropic. 
The same procedure has been carried out to elucidate 
the atomic probability distribution on the (In,Ga) sub- 
lattice. The results are shown in Fig. [l0| plotted on the 
same scale as in Fig. ||. In contrast to the As atom static 
distribution, the (In,Ga) probability distribution is much 
more isotropic and sharply peaked in space around the 
virtual crystal lattice site. 

In all compositions, the As atom distribution is highly 
anisotropic as evident in Fig. ^ with large displacements 
along (100) and (111) directions. This can be under- 
stood easily within cluster model as we discussed in Sec- 



FIG. 10. Iso-probability surface for the ensemble averaged 
(In,Ga) atom distribution. The surfaces plotted all enclose 
the volume where As atoms will be found with 68 % probabil- 
ity, (a) Ino.17Gao.s3As (b) Ino.33Gao.67 As (c) Ino.50Gao.50As 
(d) Ino.s3Gao.17As. In each case, the probability distribution 
is viewed down the [001] axis. These surfaces are plotted on 
the same scale as those in Fig. 9. 



tion IV A. The (100) displacements occur in type III 
clusters and the (111) displacements occur in type II and 
IV clusters. This also explains why, in the gallium rich 
alloy in which the three and four Ga cluster is dominant, 
the major As atom displacements are along [111], [111], 
[111], and [111] as we observed in Fig. 0(a). On the con- 
trary, in the indium rich alloy, the major displacements 
are along [111], [111], [111], and [HI], as can be clearly 
seen in Fig. |](d). 

The atomic probability distribution obtained from the 
Kirkwood model for the (In,Ga) sublattice is shown in 
Fig. [h| As we discussed, this is much more isotropic 
(though not perfectly so), and more sharply peaked than 
the As atom.. distribution. However, contrary to earlier 
predictions ,Ej and borne out quantitatively by the su- 
percell modeling, there is significant static disorder as- 
sociated with the (In,Ga) sublattice. In order to com- 
pare the magnitude of the static distortion of the (In,Ga) 
sublattice with that of the As sublattice, we calculated 
the standard deviation, a, of the As and (In,Ga) atomic 
probability distributions. This was calculated using o~i = 

x,y,z), 



(c?i(fc)) 2 , (i = x,y,z), where di refers to the 

displacement from the undistorted sublattice of atoms in 
the model supercell in x, y, and z directions, and N is 
the total number of atoms in the supercell. Table | sum- 
marizes the values of a for the As and (In,Ga) atomic 
probability distributions in the alloys. It shows that for 
all compositions the static disorder on the (In,Ga) sub- 
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TABLE I. Standard deviation of the As and (In,Ga) atom 
distributions in m^Gai-a; As alloys obtained from the Kirk- 
wood model. The numbers in parentheses are the estimated 
error on the last digit. For both As, and (In,Ga) atoms, 
a — o x = a v = <7 Z . See text for details. 



Diffuse scattering of in^Gn^ ^Ab ; 1=0,0 





x=0.17 


x=0.33 


x=0.50 


x=0.83 


a(As) (A) 
a{In,Ga) (A) 

<r(In,Ga) 
cr(As) 


0.072(1) 
0.044(1) 
0.61 


0.092(1) 
0.058(1) 
0.63 


0.097(1) 
0.060(1) 
0.62 


0.074(1) 
0.048(1) 
0.61 



lattice is around 60% of the disorder on the As sublattice. 
These static distortions give rise to a broadening of PDF 
peaks as described in Ref. 27 and evident in Fig. [5] of this 
paper. To evaluate the static contribution to the PDF 
peak broadening, ap, from the a's reported in Table || 
we used the following expression: 



(4) 



where a, b can be As, or (In,Ga). For example, for x = 0.5 
alloy, we get <j 2 D = 0.0188(4) A 2 for As-As peaks in the 
PDF, 0.0130(4) A 2 for As-(In,Ga) peaks and 0.0072(3) 
A 2 for (In,Ga)-(In,Ga) peaks. These values are in good 
agreement with the mean square static PDF peak broad- 
ening of As-As, As-(In,Ga) and (In,Ga)-(In,Ga) peaks, 
shown in Fig. 4 of Ref. 27, of 0.0187(1) A 2 , 0.0128(1) A 2 , 
and 0.0053(1) A 2 respectively. 



V. CORRELATED ATOMIC DISPLACEMENTS 

We have shown that on the average, atomic displace- 
ments of As atoms in In x Gai_ x As alloy are highly direc- 
tional. In this section, we would like to address the ques- 
tion whether these atomic displacements are correlated 
from site to site. To investigate this we have calculated 
theoretically the diffuse scattering intensity which would 
be obtained from the relaxed Kirkwood supercell model 
and compare it with the known experimental diffuse scat- 
tering. 

The Fig. [ll] shows diffuse scattering of Ina^Gao.sAs al- 
loy calculated using the DISCUS programO In this cal- 
culation the Bragg-peak intensities have been removed. 
Strong diffuse scattering is evident at the Bragg points 
in the characteristic butterfly shape pointing towards the 
origin of reciprocal space. This is the Huang scattering 
which is peaked close to Bragg-peak positions and has 
already been worked out in detail.c3 

In addition to this, clear streaks are apparent run- 
ning perpendicular to the [110] direction. The dif- 
fuse scattering calculations on (hkl) planes where I 7^ 
0, integer (Fig. [l2]), show that these diffuse streaks are 
extended along the [00Z] direction consisting of sheets 
of diffuse scattering perpendicular to the [110] direc- 
tion of reciprocal space. Diffuse scattering with exactly 
this (110) symmetry was observed in the TEM study of 
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FIG. 11. Single crystal diffuse scattering intensity ob- 
tained from the relaxed supercell model for the Ino.5Gao.5As 
alloy. The cut shown is the diffuse intensity expected in the 
(hkO) plane of reciprocal space. Bragg peaks have been re- 
moved for clarity. See text for details. 

Diffusa scattering of tag jRa^ gis ; 1=0.5 
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FIG. 12. Single crystal diffuse scattering intensity ob- 
tained from the relaxed supercell model for the Ino.5Gao.5As 
alloy. The cut shown is the diffuse intensity expected in the 
(hk0.5) plane of reciprocal space. Bragg peaks have been re- 
moved for clarity. See text for details. 



Ino.53Gao.47As.EII Careful observation of our calculated 
diffuse scattering indicates that the diffuse scattering has 



a maximum on the low-Q side of the (hkQ) planes pass- 
ing through the Bragg points, with an intensity minimum 
on the high-Q side of these planes. This is characteris- 
tic size-effect scattering obtained from correlated atomic 
displacements due to a mismatch between chemically dis- 
tinct species as recently observed ip-,a single-crystal dif- 
fuse scattering study on Si\- x Ge x ^ for example. This 
asymmetric scattering was clearly observed4n the earlier 
diffuse scattering study on Ino.53Gao.47As.El 

The single-crystal diffuse scattering intensity which is 
piled up far from the Bragg-points is giving informa- 
tion about intermediate range ordering of the atomic 
displacements. It is interesting that it is piled up in 
planes perpendicular to [110] whereas the local atomic 
displacements are predominantly along (100) and (111) 
directions. This observation underscores the complemen- 
tarity of single-crystal diffuse scattering and real-space 
measurements such as the PDF. The real-space measure- 
ments are mostly sensitive to the direction and magnitude 
of local atomic displacements and less sensitive to how the 
displacements are correlated over longer-range (though 
this information is in the data). On the other hand, the 
single crystal diffuse scattering immediately yields the 
intermediate range correlations of the displacements but 
one has to work harder to extract information about the 
size and nature of the local atomic displacements. Used 
together these two approaches, together with XAFS, can 
reveal a great deal of complementary information about 
the local structure of disordered materials. 

The single crystal diffuse scattering suggests that 
atomic displacements are most strongly correlated (i.e., 
correlated over the longest range) along [110] directions 
although the displacements themselves occur along (100) 
and (111) directions. The reason may be that the zinc- 
blende crystal is stiffest along [110] directions because 
of the elastic anisotropy in the cubic crystal. This was 
shown for the case of InAs and was used to explain why 
the 5th peak in the PDF (coming from In- As next neigh- 
bor correlations along [110] direction) was [anomalously 
sharp in both experiments and calculations.!^ If the ma- 
terial is stiffer in this direction, one would expect that 
strain fields from displacements will propagate further in 
these directions than other directions in the crystal cor- 
relating the displacements over longer range. These are 
consistent with the displacement pair correlation func- 
tion calculation by Glasca which shows that the corre- 
lation along (110) directions is larger than correlations 
along (100) and (111) and extends further. 

VI. CONCLUSIONS 

In conclusion, we have obtained high real-space resolu- 
tion PDFs of In x Gai_ x As (0 < x < 1) alloys using high 
energy synchrotron x-ray diffraction. For this purpose, 
we developed a data analysis technique adequate for high 
energy synchrotron x-ray diffraction. The PDFs show a 



clearly resolved doublet corresponding to the Ga-As and 
In-As bond lengths in the first peak of the alloys. Far- 
neighbors peaks are much broader in the alloys than that 
of the pure end members. 

We show that As atom displacements are highly direc- 
tional and can be represented as a combination of (100) 
and (111) displacements. On the contrary, the (In,Ga) 
atomic distribution is much more isotropic. The mag- 
nitude of (In,Ga) sublattice disorder is less than, but 
rather comparable (<7(i n ,Ga) ~ OS&As) to, the As sublat- 
tice disorder. Also, the single crystal diffuse scattering 
shows that atomic displacements are correlated over the 
longest range in [110] directions although the displace- 
ments themselves occur along (100) and (111) directions. 

All Oj£-|the available data,— including previous XAF-S, 
studies,!^ the present data,EZl differential PDF data,c3 
and diffuse scattering on a closely related systen£3 are 
well explained by a relaxed supercell model based on 
the Kirkwood potentialo This study also underscores 
the importance of having data from complementary tech- 
niques when studying the detailed structure of crystals 
with significant disorder. 
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